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ABSTRACT 

The method of NRQED is an adaptation of QED to bound systems. While it is fully 
equivalent to QED, it enables us to explore the QED bound states systematically as an 
expansion in the fine structure constant a and velocity (~ Za) of the bound electron. I 
describe how to construct the NRQED Hamiltonian choosing recent works on the a{Za), 
a^lZa), a{ZaY, a{Za)^, and a‘^{Za)‘^ radiative corrections to the hyperfine structure of 
the muonium ground state as examples. 


1. Introduction 


The most widely used approach to the relativistic bound state problem is that 
based on the Bethe-Salpeter equation. However, this tends to become very cumbersome 
when applied to higher-order radiative corrections, while the result reflects mostly the non- 
relativistic features of the system, which are well described by the Schrodinger equation. 
The problem is how to relate the latter to the relativistic quantum held theory. In the 
case of the electromagnetic interaction, a solution was given by the theory called non- 
relativistic quantum electrodynamics, or NRQED0. 

In this talk, I should like to describe the main features of NRQED choosing as 
an example our recent work on the radiative corrections to the hyperfine structure of 
the muonium ground state. This quantity, which is one of the most precisely measured 
quantities, is particularly suitable for a detailed examination of NRQED, being relatively 
free from the effects of hadronic interaction. The latest measurement gives @ 


Az/(exp) = 4 463 302.776 (51) kHz (11 ppb), (1) 

which includes the earlier result i and reduces its uncertainty by a factor of three. 
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2. The NRQED method 


The NRQED is a rigorous and systematic adaptation of QED to bound states. It 
is described by a Hamiltonian derived from QED by expanding it in the velocity of the 
bound electron. The equivalence of NRQED and QED is guaranteed by requiring that 
the NRQED scattering amplitude coincides with the QED scattering amplitude at some 
chosen momentum scale, e.g., at the threshold of the external on-shell particles. The 
coupling constant a and mass, renormalized on-shell, are identihed with the observed 
values and used as the input coupling constant and mass of NRQED. 

The NRQED Hamiltonian consists of all possible local interactions satisfying the 
reqnired symmetries, snch as Galileian invariance, gange invariance, parity invariance, 
time reversal invariance, hermiticity, and locality. In the absence of dynamical photons it 
is nothing bnt the Foldy-Wonthnysen-Tani (FWT) transform of the Dirac Hamiltonian. 
Since the photon is always relativistic, its Lagrangian is the same as that of QED. When 
the interaction with photons is taken into acconnt, the NRQED Hamiltonian breaks up 
into two parts: FTmain R^ntact- -^main bilinear in fermion operators as well as photon 
operators. In addition, new photon interaction terms, snch as vacnnm polarization and 
light-by-light scattering, are introdnced to represent insertion of the fermion loop. R^ntact 
represents the remainder. 

The crucial feature of NRQED is that all operators (terms of the Hamiltonian) 
are restricted to small momentum transfer, dehned by some cnt-off A, while effects of 
large momentnm transfer are represented by coefficients of these operators and terms of 
-^contact- Since this theory is meant to apply to non-relativistic systems, the cnt-off A may 
be chosen as typical mass scale of the system, e.g., the rest mass of the electron. 

For the calculation of non-recoil terms in which we can take the limit m/M 0, 

relevant to the mnoninm calculation can be written, to order a, as 
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where ip is the Pauli two-component spinor, p' and p are the ontgoing and incoming 
fermion momenta, respectively, and g = {q^,q) is the incoming photon momentnm. In the 
seagnll vertex (on the third line), ql is the incoming momentnm of the vector potential 
A. The large momentnm-transfer effect is represented by the NRQED “renormalization” 
coefficients which, to order a, are fonnd to be 
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Note that q’s do not have coefficients involving Za caused by the binding effect because 
they are determined solely by comparison of the NRQED and QED scattering amplitudes 
without referring to bound states. 

-^main aloue cauuot reproduce the amplitudes of QED. To make NRQED fully 
equivalent to QED, we must add terms of contact interaction type: 


H contact = - 

“ - ••• , ( 4 ) 


where y represents the muon (of mass M). The coefficients di are determined completely 
by the requirement that the contact terms make up the difference between the QED 
electron-muon scattering amplitude and the corresponding NRQED scattering amplitude 
generated by Tfmain lower order Tf^ntacf Several terms of TfOTntact have been deter¬ 
mined: 
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The Za factor in (^ comes from exchange of photons between the electron and the muon 
(of charge —Ze) and has nothing to do with the bound-state effect. 

Once the NRQED Hamiltonian is constructed the corrections to the energy and 
wave function can be determined by the usual bound state perturbation theory. 


3. Outline of non-recoil term evaluation 

Thus far, non-recoil terms of order a{Za), a{Za)‘^, and a‘^{Za) have been evalu¬ 
ated explicitly in NRQED. They agree with evaluation by other methods. Evaluation of 
a{Za)^ and c^iZa'f' terms within the framework of NRQED is in progress. Instead of 
discussing these calculations in detail, I shall focus on their qualitative features to illus¬ 
trate how the NRQED method works. For simplicity let us discuss here only diagrams 
with radiative photons, leaving out the vacuum-polarization contribution. 

In general, radiative corrections of order oH^iZaY' to the hyperfine splitting are of 
two types: (H) First-order perturbation of contact terms of order d^iZa)'^^^. (B) Higher- 
order perturbation of lower-order terms of Tfmain ^md 77(^ntacf Th® type (H) originates 
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from a set of QED scattering diagrams in which n + 1 photons are exchanged between 
the electron and the muon and the electron line is dressed by m virtual photon lines. 

Since operators of the NRQED Hamiltonian describe soft (low momentum-transfer) 
effects only, the type {B) scattering amplitude will reproduce the low energy behavior 
(including renormalization effects) but not the complete QED scattering amplitude. The 
difference between the QED scattering amplitude and the NRQED scattering amplitude 
is due to large momentum-transfer effect and leads to the type (H) contact term. 

The a(Za) term. The Hamiltonian consists only of terms with even parity. 

Their expectation values with respect to the Coulomb wave function are even in Za, too. 
Thus type {B) terms cannot generate terms proportional to Za. The only contribution is 
from the type (H) term, which arises from the difference of QED and NRQED scattering 
amplitudes: 

The QED amplitude is derived from the electron-muon scattering diagrams in which 

two photons are exchanged between the electron and the muon and the electron line is 
dressed by a virtual photon line. Both and have threshold singularities 

in the limit of vanishing Coulomb photon momentum, which cancel in the difference (^. 

The origin of the threshold singularity may be understood as follows: In order to 
contribute to the hyperhne splitting, one of the two exchanged photons in the bound state 
perturbation theory must be Coulomb-like while the other is transverse. The Coulomb 
photon can be absorbed in the wave function or the Green function. Such a diagram is 
reduced to lower order in Za, or, equivalently, multiplied by 1/(Zq;). In the corresponding 
scattering state, this factor manifests itself as an IR-divergent factor. Since the bound 
state has no threshold singularity all these singularities must cancel out in the end. 

The expectation value of the contact term (j^) can be evaluated in the hrst order 
bound state perturbation theory. Only the value of the Coulomb wave function at the 
origin contributes to the integral, leading to 

Ai/ = (7) 

where |0(O)p = y^/vr for the ground state. The difference between J = 1 and J = 0 spin 
states is taken. Together with the vacuum-polarization contribution this gives the hrst 
term of di in (j^). 

The a'^{Za) term. For the same reason as in the a{Za) case this correction 
comes only from the contact term of type (H). It is dehned as the difference between 
the gauge invariant set of QED scattering diagrams in which the electron and the muon 
exchange two photons and the electron line is dressed by two radiative photons and the 
corresponding NRQED scattering diagrams. 

The a{ZaY term. One contribution comes from a type (H) contact term of order 
a{ZaY, which arises from QED scattering diagrams that exchange three photons between 
the electron and the muon, dressed by one radiative photon on the electron line. These 
amplitudes have threshold singularity of up to order (m/A)^. The type {B) NRQED 
scattering diagrams which have the same leading threshold singularity are of two kinds: 
One is the hrst term of di in (|^) combined with one Coulomb potential and the other is 
the Fermi potential times the anomalous magnetic moment combined with two Coulomb 
potentials. 
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The a^{ZaY term. The relevant QED diagrams are those that exchange three 
photons between the electron and the mnon, dressed by two radiative photons on the 
electron line. The corresponding NRQED diagrams have two radiative photons, too. 
However, these diagrams contribnte only to order a‘^{Za)^, which is too small by a factor 
(Za)^. Thns the radiative correction of order a‘^{Za)^ comes from the QED diagrams 
only. The contact term of NRQED derived from these QED diagrams does not have the 
ln(m/ 7 ) dependence. Thns logarithmic terms of order can come only from the 

NRQED operators of contribnting to the a(Za)^ term combined with an additional 
factor of a fonnd in the “renormalization” effect. This leads to the resnlt reported in (|^. 

The a(Za)^ term. The status of the theory of this term will be discussed in 
detail by Dr. Nio in her talk presented at this Workshop. 


4. Current status of theory of muonium hfs 

As is well known, the bulk of the hyperhne splitting is given by the Fermi formula 
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= 4 453 839.405 (518) kHz, 


( 8 ) 


where Roo is the Rydberg constant for inhnite nuclear mass, and m and M are the dectron 
and muon masses, respectively. The values used above for a, Roo, and M/m are 


= 137.035 999 58 (52) (3.8 ppb), 

Roo = 10 973 731.568 639 (91) m-\ 
M 

— = 206.768 273 (24) (117 ppb). 
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Various theoretical corrections to Ep have been calculated over the last 50 years. 
They may be classihed into non-recoil and recoil terms, both including radiative and 
binding effects. In addition there are contributions of hadronic vacuum-polarization and 
weak interaction effect. Thus one may write 

Az/(theory) = An(non-recoil) -|- Az/(recoil) -|- Ai/(hadron) -|- Az/(weak). (10) 

A purely radiative non-recoil term of order a{Za) and an approximate value of the 
a{ZaY term have been known for some time i. Recently radiative corrections of orders 
a(Za)^ and a‘^{Za) have been evaluated with high precision. Including these results and 
the Fermi term Ef itself, the non-recoil term of the muonium hyperfine structure can be 

written asBSi SB 0 
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Here and are the anomalous magnetic moments of the electron and muon, respec¬ 
tively, and b is the Breit term 
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The appearance of the factor (1 -|- a^) in ([Tl| ) is in accord with the dehnition (B) of Ep. 
The coefficient 16.9042 of the a{ZaY term is different from 15.39 given in H. This is 
because the latter contains, besides the a{ZaY term, terms of order a(Za)^ and higher. 
More recent evaluation of this term without expansion in Za is given in li3. For a detailed 
discussion see Ii3. 

In order to match the precision of the new measurement H, however, it is necessary 
to improve the theoretical prediction of terms of order This is because, when enhanced 
by the factor ln^(m/ 7 ) or ln(m/ 7 ) where m /7 = 1 /Za, corrections of order a{ZaY and 
a‘^{ZaY contribute to the muonium hyperhne structure as much as the a^{Za) term. 
(There is no In(Za) enhancement for a^{Za) and {ZaY terms.) 

The leading log contribution of the a{ZaY term is known exactly H 0 HI 0. The 
non-log term has also been estimated 0. Altogether we have 
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where the hrst term in the brackets is the contribution from the radiative photon and the 
second is due to the vacuum polarization. The uncertainty is an estimate based on the 
non-log part of the a(Za)^ term given in 0. We are trying to improve this uncertainty. 

For the a^fZaY term the leading log part together with an estimate of non-log 
part gives the value 0 ’ 10 
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where the factor m/M originates from the reduced mass effect. An error in B due to an 
incorrect treatment of low energy contribution is rectihed in ([T^). 

The pure recoil and order a radiative recoil corrections, except for the last term of 
order {m/MY obtained recently 0 , have been known for some time i’ 0 0 


Az/(recoil) = Ei 
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where 7 = Zarrir, rrir = mM/{m + M). The In^ and In^ parts of the a^iZot) term were 
evaluated by Hides et al. lil. 

Further rehnement of A(recoil) comes from evaluation of {Za)^{m/M) terms. One 
of the leading log term of this contribution arises from the recoil hfs potential (the hrst 
line of (|T^) and the 1-photon exchange part of the Coulomb Green function: 

Az/«((Za)3(m/M)) = 2 < YmsG^^Vk >+2 < YmsC.^Vd > 


3{Za)^ m 


TT M 

= -0.210 (43) kHz, 


(-) 

\m / 


Erlni — In 


7 


(16) 


where Yk and Yp are the kinetic energy and Darwin terms, respectively. The uncer¬ 
tainty is estimated assuming that it is a factor ln{Za) smaller than the leading term. 

Another log term arises from the Salpeter recoil correction to the Lamb shift i^. 
This contribution involves both relativistic and non-relativistic regions. Unfortunately, 
the latter was handled incorrectly in i. Correcting this error we hnd 
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where one factor of Za is radiative correction due to exchange of transverse photon 
between the electron and the muon and Cs is given by 

Cs = -- + -(21n2 -h 3) - - 


—77-77:7 flii(l + m/M) — im/MY ln(l -|- M/m)) . (18) 
[m/M)^ ^ > 


The hadronic vacuum polarization contributes hi 

Ai/(hadron) = ^^^^^(3.75 ± 0.24)^^ 
tt" m" 

= 0.250 (16) kHz, 
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where m7r is the charged pion mass. 
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Finally, there is a contribution due to the exchange 

A / 1 \ ^ 3y/2mM 

Ai/(weak) = —Gp - Gp 

San 

~ —0.065 kHz, 


where Gp is the Fermi weak coupling constant. 

Collecting all contributions we obtain the best current estimate 

Az/(theory) = 4 463 302.649 (517) (34) (< 100)kHz , 


( 20 ) 


( 21 ) 


where the hrst and second errors come from those of M/m and a in (|^) and the last is 
an estimate of the theoretical uncertainty from the yet-to-be calculated terms. This is in 
excellent agreement with the measurement (|ID. 


5. Concluding remarks 


Current uncertainties in (|^, (P!4|), (PIB]), and (P!7|) are 0.008 kHz, 0.024 kHz, 0.043 
kHz, and 0.059 kHz, respectively. For further improvement it is particularly important to 
evaluate the 0{a^{m/M)) recoil terms. 

Study of positronium in the framework of NRQED requires additional terms to the 
Hamiltonian. This is discussed by P. Labelle in this Workshop. Adaptation of NRQED 
method to the chiral perturbation theory is of great interest at present in view of the hne 
experimental works on the bound system. 

This work is supported in part by the U. S. National Science Foundation. 
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